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I. INTRODUCTION 



The surprising discovery of new states decaying to 
has renewed interest in the field of 
charmonium spectroscopy, since not all the new reso- 
nances are easy to accommodate in the quark model. 
Specifically, the BABAR experiment discovered a broad 
state, A (4260), decaying to J/ipn + Tr~, in the initial- 
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state-radiation (ISR) reaction e+e~ — > 7/si?A"(4260). 
Its quantum numbers J PC = 1 are inferred from 
the single virtual-photon production. Enhancements in 
the tl)(2S)n + TT- mass distribution at 4.36 GeV/c 2 0, 0] 
and 4.66 GeV/c 2 Q have been observed for the reac- 
tion e + e~ — > r )isR^p{2S)'K + 'K^ . Charmonium states 
at these masses would be expected || |t| to decay 
predominantly to DD, D*D, or D*D* . It is pecu- 
liar that the decay rate to the hidden charm final 
state J/-07r + 7r~ is much larger for the X(4260) than 
for the higher-mass charmonium states (radial excita- 
tions) 0j- Many theoretical interpretations for the 
A (4260) have been proposed, including unconventional 
scenarios: quark-antiquark gluon hybrids @, baryo- 
nium [Tol | . tetraquarks [111) and hadronic molecules [l^ . 
If the X(4260) were a diquark-antidiquark state [cs][cs], 
as proposed by L. Maiani et al. [ill ], this state would pre- 
dominantly decay to D+Dj. For a discussion and a list 
of references see, for example, Ref. 
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In this paper, we present a study of the ISR production 
of D+D~, D*+D-, and D*+ D*~ Q final states, and 
search for evidence of charmonium states and resonant 
structures. This follows earlier BABAR measurements of 
the cross section of DD [15[ and of D*D and D*D* pro- 
duction [r| and studies of these final states [13, 13 
the Belle Collaboration. Recently the CLEO Collabora- 
tion [19| studied e+e~ annihilation to D+D~, D* + Dj, 
and D* s + D* s ~~ final states at center of mass energies from 
threshold up to 4.3 GeV GeV/c 2 . In the present analysis 
we extend these measurements up to 6.2 GeV/c 2 . 

This paper is organized as follows. A short description 
of the BABAR experiment is given in Section II, and data 
selection is described in Section III. In Sections IV, V, 
and VI we present studies of the DfD~, D* s + Dj , and 
D* + D*~ final states, respectively. Fits to the three fi- 
nal states are described in Section VII, and in Section 
VIII we present limits on the decay of the V(4260) to 
summary and conclusions are found in 

Section IX. 



II. THE BABAR EXPERIMENT 

This analysis is based on a data sample of 525 fb _1 
recorded mostly at the T(4S I ) resonance and 40 MeV 
below the resonance by the BABAR detector at the PEP- 
II asymmetric-energy e + e~ storage rings. The sample 
includes also 15.9 fb _1 and 31.2 fb _1 collected at the 
T(2S*) and T(35) respectively, and 4.4 fb _1 above the 
T(4S') resonances. The BABAR detector is described in 
detail elsewhere [20| . We mention here only the compo- 
nents of the detector that are used in the present analysis. 
Charged particles are detected and their momenta mea- 
sured with a combination of a cylindrical drift chamber 
(DCH) and a silicon vertex tracker (SVT), both oper- 
ating within a 1.5 T magnetic field of a superconducting 
solenoid. Information from a ring-imaging Cherenkov de- 
tector is combined with specific ionization measurements 
from the SVT and DCH to identify charged kaon and 
pion candidates. The efficiency for kaon identification 
is 90% while the rate for a kaon being misidentified as a 
pion is 2%. Photon energies are measured with a CsI(Tl) 
electromagnetic calorimeter (EMC). 

III. DATA SELECTION 

For each candidate event, we first reconstruct a D+D~ 
pair. While one of the Df is required to decay to K + 
K~ 7r + , we include three different decay channels for the 
second Dj (see Table fl]) . D* + decays are reconstructed 
via their decay D* + -> D+7. 

For all final states, events are retained if the number 
of well-measured charged tracks having a transverse mo- 
mentum greater than 0.1 GeV/c is exactly equal to the 
total number of charged daughter particles. EMC clus- 
ters with a minimum energy of 30 MeV that are not as- 



TABLE I: Reconstructed decay channels for the two 
mesons in each event. 

Channel First D s decay Second D s decay 

(1) K+ K~ 7T+ K+ K~ tt- 

(2) K+ K- 7T+ K+ K' 7r" 7T° 

(3) K+ FT 7T+ K% K~ 



sociated with a charged track are identified as photons. 
Candidates for the decay tt° — > 77 are kinematically con- 
strained to the 7T° mass. For K® —> 7r + 7r~ candidates 
we apply vertex and mass constraints. The tracks cor- 
responding to the charged daughters of each D+ candi- 
date are constrained to come from a common vertex. Re- 
constructed Dj candidates with a fit probability greater 
than 0.1% are retained. Each DfD~ pair is refit to a 
common vertex with the constraint that the pair origi- 
nates from the e + e~ interaction region. Only candidates 
with a x 2 fit probability greater than 0.1% are retained. 
For each event we consider all combinations. 

ISR Monte Carlo (MC) [U events for each final state 
are fully simulated using the GEANT4 detector simula- 
tion package [22j], and they are processed through the 
same reconstruction and analysis chain as the data. 

We select D+ and D* + candidates using the re- 
constructed mass and the mass difference, which 
for D+ K + K~tt + is defined as Am(Dfj) = 

m(K + K~n + ^) — m(K + K~ir + ). The Df parameters 
are obtained by fitting the relevant mass spectra using 
a polynomial for the background and a single Gaussian 
for the signal. For D* + , we use the PDG [8| mass and a 
Gaussian width a — 6 MeV/c 2 obtained by MC simula- 
tions. Events are selected within ±2. Oct from the fitted 
(*)+ 

central values. The D s candidate three-momentum is 
determined from the summed three-momenta of its decay 
particles. The nominal mass || is used to compute 

the energy component of its four-momentum. 

The ISR photon is preferentially emitted at small an- 
gles with respect to the beam axis, and it escapes de- 
tection in the majority of ISR events. Consequently, 
the ISR photon is treated as a missing particle. We de- 
fine the squared mass M r 2 ec recoiling against the D+D~ , 
D* + D~, and D* + D*~ systems using the four-momenta 

of the beam particles p e ± and of the reconstructed Ds 

P D M±- 

M rec = (Pe- + Pe+ ~ P D M+ ~-Pd<*'-) 2 ' W 

This quantity should peak near zero for both ISR events 
and for exclusive production of e + e~ — > D y s ' D y s ' . 
For exclusive production, the D s D s mass distri- 
bution peaks at the kinematic limit. We reject exclu- 
sive events by requiring the D s *^ D s *^~ mass to be be- 
low 6.2 GeV/c 2 and select ISR candidates by requiring 
|M 2 ec | < 0.8 GeV 2 /c 4 . 
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FIG. 1: Distributions of M r 2 ec for the (a) D+D~, (b) D* s + D~ , and (c) D* + D*~ final states. 



We allow additional tt° and photon candidates due to 
radiative or background photons. This introduces multi- 
ple candidates in the reconstruction of the different chan- 
nels. For channel (l)-(2) ambiguities, each [K + K~ tt + 
(tt )} combination is considered as a candidate for both 
channels. For all channels, each [Df, Dj , (7), (7)] com- 
bination is considered as a candidate for DfD~, D* s + Dj , 
and D* s +D* s -. 

To discriminate among the different Df channels and 
D s ' D s ' final states, and to separate signal from back- 
ground, we make use of a likelihood ratio test: 



N 



(2) 



where N is the number of discriminating variables, while 
PDFf and PDFf are normalized distributions describ- 
ing signal and background, respectively. Signal PDFf 
are obtained from MC simulations. Background PDFf 
are obtained from the data. Since the ISR signal is very 
small compared to the entire data set of candidates (< 0.1 
%), we use the data as the background model by relax- 
ing all the selection criteria, except m{D s *^ Jr D s *^ 1 ) < 6.2 
GeV/c 2 . 

The discriminating variables used in the likelihood ra- 
tio are the following. 

• The number of additional 7r° candidates in the 
event. For decay channel (2) this number is com- 
puted after removing the ir° from Df decay. This 
distribution is expected to peak at zero for signal 
events. 

• The residual energy in the calorimeter, which is 
computed after removing any ISR photon candi- 
date, identified by a center-of-mass energy greater 



than 2.0 GeV/c 2 . For D*+D~ and D*+D*- final 
states, the 7 from D* s + decays is excluded from the 
residual energy calculation. This distribution is ex- 
pected to peak at zero for signal events. 



• The distribution of cos#*, where 9* is the polar 
angle of the D s D s system in the center-of- 
mass frame which peaks at ±1 for ISR events. 



• The momentum distribution of the 7r° from the Df 
for decay channel (2). 



• The 7 energy distribution from D* + for 2D* 
and D* + D*~ final states. 



D: 



For each Df decay channel (l)-(3) and for each D+Dj , 
D* + Dj, and D* + D*~ final state, we produce a likeli- 
hood ratio test according to Eq. (2) and apply empirically 
determined cuts on L in order to reduce the background 
and minimize the signal loss. 



IV. STUDY OF THE DJ D~ FINAL STATE 

MC studies demonstrate that the main background to 
the D+Dj final state is from D* + Dj events, which have 
a larger cross section. Therefore, we eliminate DfD~ 
candidates if they are also identified as D* + Dj candi- 
dates. MC simulations show that this veto rejects about 
14% of the true D+D~ final states and that the residual 
background is consistent with D* + D~ feedthrough. 

FigureQJa) shows the M r 2 c distribution for the selected 
DfDj candidates, summed over the the D+ decay chan- 
nels (l)-(3). The peak centered at zero is evidence for 
the ISR process. To determine the number of signal and 
background events, we perform a \ 2 fit to the -M 2 CC distri- 
bution. The background is approximated by a 2 nd order 
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TABLE II: Number of signal ISR candidates and purities 
for the different final states calculated in the range |M 2 CC | < 
0.8 GeV 2 /c 4 . 



Final state 


Signal+ Background 


Purity (%) 


D+D- 


8f 


65.4 ± 5.3 


D*+D7 


286 


67.f ± 2.8 


D*+D*- 


f05 


54.3 ± 4.9 



polynomial. The signal lineshape is taken from DfD s 
MC simulations. The resulting yield and the fitted pu- 
rity P, defined as P = iV signal / (iV signal + ^background) are 
summarized in Table [TTJ 

The D+D~ mass spectrum, presented in Fig. [2ja), 
shows a threshold enhancement at the position of the 
-0(4040) and a small enhancement around 4.26 GeV/c 2 . 

We make use of the Gaussian functions to describe the 
presence of the peaking backgrounds. The D+D~ back- 
ground, taken from M r 2 ec sideband events (1.5 < |M 2 ec | < 
3.5 GeV 2 /c 4 ), is fitted to a sum of a Gaussian function 
and a 3 rd order polynomial. The fitted DfD~ mass spec- 
trum for these events, normalized to the background esti- 
mated from the fit to the M r 2 c distribution, is presented 
as the shaded distribution in Fig.[2ja). 

The D+D~ reconstruction efficiency and the mass res- 
olution for each channel have been studied in the mass 
region between 4.25 and 6.25 GeV/c 2 . The D+D~ mass 
resolution is similar for decay channels (1) and (3) and 
slightly worse for decay channel (2) (by w 1 MeV/c 2 ). It 
increases with DfDj mass from 3.5 to 5.5 MeV/c 2 in 
the mass region of the t/j resonances (< 5 GeV/c 2 ). The 
mass-dependent reconstruction efficiency for the DfD~ 
decay channel % (i = 1,3), ei(m D + D -), evaluated at five 
different mass values, is parameterized in terms of a 2 nd 
order polynomial, and scaled to account for the prod- 
uct branching fractions for each channel, Bi Q, given in 
Table HI 



20 
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FIG. 2: The observed (a) D+ D~ , (b) D* 3 + D~, and (c) 
Ds + D$~ mass spectra. The shaded areas show the back- 
ground derived from fits to the M r 2 ec sidebands. The dashed 
lines indicate the sum of this background and the coherent 
background. The solid lines are the results from the fit de- 
scribed in Section VII. 



i(m 



dtd: 



(3) 



These values are weighted by Ni(m D + D -), the number 
of DfD~ candidates in decay channel i, to compute the 
average efficiency as a function of m D + D - , 



\ m DiD-) 



(4) 



ef (m 



The e function for DfD~ is shown in Fig. [3ja). The 
three DfD~ decay channels, after correcting for effi- 
ciency and branching fractions, have yields that are con- 
sistent within statistical errors. 



The DfD s cross section is computed using 

dN/dm, 



1 djd: 



DID. 



)dC/dm 



dtd: 



(5) 

where dN/dm D + D - is the background-subtracted yield. 
The differential luminosity is computed as [23| 



dC 



2m 



d+d^__o_ 

1TX 



(ln(s/m 2 ) - l)(2-2a; + a; 2 ), 

(6) 



where s is the square of the e + e center-of-mass energy, 
a is the fine-structure constant, x = 1 — m 2 , „_ Is, m e 

DJD S 1 ' 

is the electron mass, and C is the integrated luminosity 
of 525 fb" 1 . The cross section for DfDj is shown in 
Fig.[UJa). This result can be compared with QCD calcu- 
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FIG. 3: Weighted efficiencies e B for (a) D+DJ, (b) D*+ D~ , and (c) D* + D* a ~. 



TABLE III: Systematic uncertainties (in %) for the evaluation 
of the Df Dj , Dg + Dj , and Dl + Dl~ cross sections. 



Source 


D+D- 


nr d- 


D* s + D*~ 


Background subtraction 


18.0 


4.2 


4.9 


Branching fractions 


10.0 


10.0 


10.0 


Particle identification 


5.0 


5.0 


5.0 


Tracking efficiency 


1.4 


1.4 


1.4 


7T°'s and 7 


1.1 


2.9 


4.7 


Likelihood selection 


8.7 


4.0 




Total 


23 


13 


13 



purity are summarized in Table [TTJ 

The D* s + Dj mass spectrum and background are 
shown in Fig. HJb) and is dominated by the ■0(4160) res- 
onance. The D* + D~ mass resolution is similar for the 
three decay channels and increases with D* + Dj mass 
from 7 to 8 MeV/c 2 in the mass region of the ip reso- 
nances. The weighted efficiency e B is shown in Fig.[3]Jb). 
The D* + D~ cross section is calculated using the method 
described in Sec. IV for DfD~. The result is shown 
in Fig. 0Jb). The overall systematic error for the cross 
section is 13% and is dominated by the uncertainties in 
the branching fractions [8j (see Table Mil) . 



lations ref. [24|, which predict a vanishing cross section 
near 5 GeV/c 2 . 

The list of systematic uncertainties for the DfD~ cross 
section is summarized in Tabic Mil and it is evaluated to 
be 23%. It includes contributions from particle identifica- 
tion, tracking, photon and 7r° reconstruction efficiencies, 
background estimates, branching fractions, and the cri- 
teria to select the final state. All contributions are added 
in quadrature. The Dj D~ systematic error is dominated 
by the uncertainty in the veto of the D* + Dj events. 



V. STUDY OF THE D* + D~ FINAL STATE 

A similar analysis is carried out for the D* + D~. Fig- 
ure Elja) shows the Am(Df~f) distributions for D* + Dj 
candidates passing the ISR requirements described in 
Sect. III. We also require the presence of a reconstructed 

The M 2 CC distribution for D* + D s candidates is shown 
in Fig. HJb) where a clear signal of ISR production is 
observed. The number of ISR candidates and sample 



VI. STUDY OF THE D* + D*~ FINAL STATE 

For the selection of D* + D*~ candidates, we do not 
make use of the likelihood test described in the previ- 
ous sections because no improvement for the signal to 
background ratio is obtained. Instead, we require the 
two photon invariant mass 771(77) to lie outside the ir 
window. 

The ambiguity in the 7's assignment to D* + or D'*~ is 
resolved by choosing the Df 7--D7 7 combinations with 

both ± masses closest the expected value. Fig. [TJc) 
shows the resulting M r 2 cc distribution which shows clear 
evidence for the signal final state produced in interactions 
with ISR. For the selected the ISR signal candidates, we 
show the Am(Df"f) distribution (two combinations per 
event) in Fig.[5]Jb). The resulting event yield and purity 
are summarized in Table [TTJ 

The mass spectrum and background are 

shown in Fig. HJc). Due to the presence of structures, 
the background in this case is fitted using a 3 rd order 
polynomial and two Gaussians. Monte Carlo studies in- 
dicate that an important part of this background is due 
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FIG. 4: Cross section for e+e~ -> (a) D+D~, (b) D* s + D^, 
and (c) Dg + Ds~. The error bars correspond to statistical 
errors only. 



to the D^ + Dj final state plus a random background 7. 
The D* + D*~ mass resolution is similar for the three de- 
cay channels. It increases with D* + Z?*~ mass from 9 to 
11 MeV/c 2 in the mass region of the ip resonances. The 
weighted efficiency e B is shown in Fig.[3fc). The D* + D*~ 
cross section shown in Fig. |4jc) is calculated using the 
same method used to compute the DfD~ cross section. 
The overall uncertainty for the cross section is 13% and is 
dominated by the uncertainties on the D s branching 
bractions (see Table UTTj). 



The D+Dj, D* + Dj, and D*+D*~ cross sections, 
where they overlap, are in good agreement with 
CLEO [ly| measurements. 
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FIG. 5: Am distributions for D* + candidates after applying 



the \M? CC \ < 0.8 GeV7c 4 and m{D { * )+ D { * ] ~) < 6.2 GeV/c 2 
selections, for the (a) D* + D~ , and (b) D* + D*~ samples. The 
shaded regions indicate the ranges used to select the D| can- 
didates. 



VII. FIT TO THE MASS SPECTRA 

Unbinned maximum likelihood fits are performed sep- 
arately to the D+D~, D* + Dj, and D* s + D* s ~ mass spec- 
tra. The likelihood function used is 



fe B (m)\P{m) + c 1 W 1 {m)e i ' t ' 1 + 



+ c n W n (m)e^\ 2 
+fl(m)(l-/),(7) 



where m is the mass, c$ and 4>i are & ee 

parameters, Wjirn) are P-wave relativistic Breit-Wigner 
distributions [a], P(m) represents the nonresonant con- 
tribution, B{m) is the background described in Sect. IV, 
e B (m) is the weighted efficiency, and / is the signal frac- 
tion fixed to the values obtained fitting the M r 2 ec distri- 
butions. In this way we allow interference between the 
resonances and the nonresonant contribution P(m). The 
shape of the nonresonant contribution P(m) is unknown; 
we therefore parametrize it in a simple way as 



P(m) = C(m)(a + bm), 



(8) 



(*)+r>(*)" 



where C(m) is the phase space function for D s D 
and a and b are free parameters. The size of the nonres- 
onant production is determined by the fit. 

The mass and width of the ^(4040), ^(4160), ii(4415) 
and X(4260) are fixed to the values reported in [8|. Res- 
olution effects can be ignored since the widths of the 
resonances are much larger than the experimental reso- 
lution. 

The three D+D~, D* + D~, and D*+D*~ likelihood 
functions are computed with different thresholds, efficien- 
cies, purities, backgrounds, and numbers of contributing 
resonances appropriate for each final state. The results 
of this fits are compared to the data in Fig. [21 both the 
total fitted yield as well as the coherent nonresonant con- 
tribution, |P(m)| 2 , ignoring any interference effects. 
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TABLE IV: D+D~, D* + D~ , and D* a + D* s ~ fit fractions (in 
%). Errors are statistical only. 



and 



Resonance 




Fraction 






DfD- 




D*+D*- 


P(m) 


11 ± 5 


27 ± 5 


71 ± 20 


0(4040) 


62 ± 21 






0(4160) 


23 ± 26 


53 ± 8 




0(4415) 


6 ± 11 


4 ± 2 


5 ± 12 


A(4260) 


0.5 ± 3.0 


18 ± 24 


11 ± 16 


Sum 


103± 36 


102 ± 26 


87 ± 28 



The fraction for each resonant contribution i is defined 
by the following expression: 



fi = 



N 2 jmjmtfdm 
^kCjclJW.MW^Mdm- 



(9) 



The fractions fi do not necessarily add up to 1 because 
of interference between amplitudes. The error for each 
fraction has been evaluated by propagating the full co- 
variance matrix obtained by the fit. The resulting fit 
fractions are given in Table HVl The D+Dj cross section 
is dominated by the ■0(4040) resonance, and the D* + D~ 
cross section by the 0(4160) resonance. The -D* + -D*~ 
cross section shows little resonance production. The fits 
to the D+D~, L»* + L»7, and D* s + D* s ~ mass spectra in- 
clude the X(4260) resonance, which is allowed to inter- 
fere with all the other terms. In all cases, the AT (4260) 
fraction is consistent with zero. We note that the weak 
enhancement around 4.26 GeV/c 2 in the D+D~ mass 
spectrum is described by the fit in terms of interference 
between the 0(4040) and 0(4160) resonances. 



VIII. LIMITS ON A(4260) 

The A" (4260) yields are used to compute the cross sec- 
tion times branching fraction, to be compared with a 
BABAR measurement of the J/iJjTr + ir~ final state 0]. The 
fractions from the fits reported in Table UVl are converted 
to yields which are divided by the mass dependent e B ef- 
ficiency and the integrated luminosity. Systematic errors 
due to the mass and the width of the 0(4040), -0(4160), 
-0(4415), and AT(4260) resonances are evaluated by vary- 
ing the masses and widths by their uncertainty in the 
fit. The size of the background contributions is varied 
within the statistical error, and the meson radii in the 
Breit-Wigner terms [25j are varied between and 2.5 
GeV -1 . Statistical and systematic errors are added in 
quadrature. We obtain 



g(X(4260) ->• D+DJ) 
£(AT(4260) -4 J/V-vr+vr-) 

£(X(4260) -> D*+D~) 
£(AT(4260) J/07T+7T-; 



< 0.7, 



< 44, 



(10) 



(11) 



i3(AT(4260) -> D*+D 



£(AT(4260) -> J/07T+7T- 

at the 95% confidence level. 



< 30, 



(12) 



IX. TOTAL CROSS SECTION AND 
CONCLUSION 

The sum of the e+e - -)• D+D~, e+e~ -)• D*+D~, and 
e+e~ — > D* + .D*~ cross sections is shown in Fig. [5] the 
arrows indicate the position of the different resonances 
and the A"(4260). At the AT(4260) mass, there is a lo- 
cal minimum, similar to the measured cross section for 
hadron production in e + e~ annihilation Q. 
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FIG. 6: Sum of e+e~ -)■ D+D~, e+e" -> D* + D~, and 
e + e~ — > Dg + Dg~ cross sections. Errors are statistical only. 
The arrows indicate the positions of the different resonances 
and the A (4260). 



In conclusion, we have studied the exclusive ISR pro- 
duction of the DfDj , D* + D~ , and D* s + D* s ~ final states. 
The mass spectra show production of the J PC = 1" 
states, 0(4040), 0(4160) and a weak indication for a 
smaller enhancement near 4.3 GeV. From fits to the mass 
spectra for the three different final states we have deter- 
mined contributions by different cc resonances. 

Upper limits on AT (4260) decays to these final states 
relative to J/ipir + ir~ are computed. If the A'(4260) is 
a 1 charmonium state, it should decay predominantly 
to open charm. Within the present limited data sample 
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size, no evidence is found for A (4260) decays to DfD s , 
D*+D-, and D*+D*-. If the X(4260) were a tetraquark 
state, it would decay predominantly to D+D~ [ill ]. 
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